Extended DFT+U+V method with on-site and inter-site electronic interactions 
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In this article we introduce a generalization of the popular DFT+U method based on the extended 
Hubbard model that includes on-site and inter-site electronic interactions. The novel corrective 
Hamiltonian is designed to study systems for which electrons are not completely localized on atomic 
states (according to the general scheme of Mott localization) and hybridization between orbitals from 
different sites plays an important role. The application of the extended functional to archetypal Mott 
- charge-transfer (NiO) and covalently bonded insulators (Si and GaAs) demonstrates its accuracy 
and versatility and the possibility to obtain a unifying and equally accurate description for a broad 
range of very diverse systems. 

PACS numbers: 71.10.Fd, 71. 15. Mb, 71.20.-b 



I. INTRODUCTION 

Due to its moderate computational cost, the simplic- 
ity of its formulation, and the ability to capture the ef- 
fects of static correlation, the DFT+U method [J, [|J (by 
DFT+U we generically indicate Hubbard-corrected Local 
Density Approximation, LDA+U, or Generalized Gradi- 
ent Approximation, GGA+U) has become quite popu- 
lar in recent years to study systems with strongly local- 
ized (typically d or /) and thus, correlated, valence elec- 
trons. DFT+U consists of a simple corrective Hamilto- 
nian, added to the (approximate) DFT energy functional, 
that is shaped on the Hubbard model [!, 0, [H, 0, 0, Q • 
In its single-band formulation the Hubbard Hamiltonian 
can be written as follows: 
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and 



7ii )Cr are electronic creation, annihilation and number op- 
erators for electrons of spin a on site i. The Hubbard 
model is normally used to describe the behavior of sys- 
tems in the limit of strong electronic localization. In 
these conditions the motion of electrons is described by a 
"hopping" process from one atomic site to its neighbors 
(first term of Eq. [1]) whose amplitude t is proportional 
to the bandwidth of the system and thus represents the 
single-particle term of the total energy. The Coulomb 
repulsion, instead, is only accounted for between elec- 
trons on the same atom through a term proportional to 
the product of the occupation numbers of atomic states 
on the same site. The effective strength of the "on-site" 
Coulomb repulsion is U. The explicit account for on-site 
electronic repulsions is ideally suited to study Mott insu- 
lators. In fact, in these systems the insulating character 
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of the ground state stems from the dominance of short- 
range Coulomb interactions (the energy cost of double 
occupancy of the same site) over single-particle terms of 
the energy (generally minimized by electronic dereal- 
ization) [9| that leads to the localization of electrons on 
atomic orbitals. 

Commonly used approximations to the exact DFT en- 
ergy functional generally describe electron-electron inter- 
actions using a mean-field-like approach (i.e., through 
functionals of the electronic charge density) which is 
not ideally suited to represent the Coulomb repulsion 
between localized electrons nor their "correlated" mo- 
tion. Furthermore, the presence of residual electronic 
self-interactions in the approximate DFT functionals, en- 
hances the tendency towards derealization. As a result, 
most approximate functionals fail to capture the behav- 
ior of systems whose electrons are strongly localized. In 
spite of its simplicity, the Hubbard-rooted DFT+U ap- 
proach has been very effective in alleviating these well 
known difficulties and has been successfully applied to 
the study of a large number of quite diverse systems. 
In recent years it has also been used as the starting ap- 
proximation for more sophisticated computational meth- 
ods as, for example, the Dynamical Mean-Field Theory 
(DFT+DMFT) [f^[Ti| and, more recently, the GW ap- 
proximation [13 . 1 1 31 ] . 

The "on-site" Hubbard Hamiltonian contains all the 
essential ingredients to capture the physics of Mott lo- 
calization in strongly correlated systems and, in fact, 
a good description of these systems is usually obtained 
with DFT+U. A generalized ("extended") formulation 
of the Hubbard model, with Coulomb interaction terms 
between electrons on neighbor sites, has been considered 
since the early days of the Hubbard model @, Q • How- 
ever, to the best of our knowledge, it has never been 
implemented in any DFT-based functional, nor used in 
ab-initio calculations. A simplified formulation of the 
original extended Hamiltonian (with only charge interac- 
tions between couples of sites) is given in the following 
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expression: 

H Hu b = t ^ (q |0 -Cj,<r + fe-c) + U ^ ^,t n U 

(i.j>><» « 

+ Vj2^ n i (2) 

where V represents the strength of the interaction be- 
tween electrons on neighbor sites i and j. 

The use of the extended model has been stimulated 
in the last decades by the discovery of high T c supercon- 
ductors and the intense research activity focusing around 
them. Whether the inter-site coupling V has a determi- 
nant role in inducing superconductivity is still matter 
of debate. The "resonating valence bond" model [l4| 
predicts a superconducting state (at least within mean- 
field theor y) f or a doped Mott insulator with only on-site 
couplings [15| . However, several numerical studies sug- 
gest that the inter-site interaction plays indeed an im- 
portant role [la . |17| and superconductivity is predicted 
in a regime with repulsive on-site (U > 0) and attrac- 
tive inter-site (V < 0) couplings pj| GJ, [M EJ. Sev- 
eral studies have also demonstrated that in the "normal" 
(non-superconducting) state of superconductors and, in 
general, in correlated materials, the relative strength of 
U and V controls many properties of the ground state 
as, for example, the occurrence ofpossible phase sep- 
arations [12], the magnetic order [23|, [24j, the onset of 
charge-density and spin-density wave regimes [25| ■ More 
recently, the extended Hubbard model has been used to 
study the conduction and the structural properties of 
polymers and carbon nano-structures and the interplay 
between U and V was shown to control, for example, the 
dimerization of graphene nanoribbons [26| . The impor- 
tance of a more accurate account of inter-site couplings 
and of the detailed balance with on-site interactions has 
been also recognized in a theoretical study of Fe impuri- 
ties deposited on various alkali metal films (27[. In this 
case, however, the hybridization between the impurity 
d states and the s and p orbitals of the alkali metal 
substrate was described using a one-body term in the 
Hamiltonian (effectively corresponding to a generalized 
hopping process between the localized and the delocal- 
ized states) rather than a two-body inter-site electronic 
interaction. 

In this work we propose to introduce the extended 
Hubbard model in DFT-based calculations through a 
generalization of the corrective functional DFT+U is 
based on to include on-site and inter-site electronic in- 
teractions. The main aim of the novel formulation is to 
improve the accuracy of the DFT+U scheme and to make 
it quantitatively predictive for systems where correlation 
is not strong enough to induce a complete Mott localiza- 
tion of electrons on atomic states, or for which, in general, 
the hybridization between orbitals belonging to different 
atoms plays an important role in determining the prop- 
erties of the ground state. The novel functional is tested 
on quite diverse bulk solids as NiO, Si and GaAs. In 



NiO, a prototype strongly correlated material, d states 
are quite well localized around Ni atoms. Nevertheless 
their hybridization with O p states plays a quite signif- 
icant role and is one of the main factors to determine 
the charge-transfer insulating character of the material. 
Si and GaAs are, on the opposite extreme, typical band 
insulators. The hybridization between s and p orbitals 
belonging to the same and to neighbor sites plays a dom- 
inant role in determining the electronic structure of these 
materials as it leads to the formation of bonds in between 
atoms (with tetrahedral symmetry) and to the conse- 
quent onset of the semiconducting character of these sys- 
tems due to the splitting between bonding (valence) and 
anti-bonding (conduction) states. The choice of these 
systems is not casual: a computational scheme that is 
able to describe accurately strongly localized as well as 
"strongly hybridized" systems is likely to be successful 
for intermediate situations as well that are, by far, the 
most difficult to deal with. For each of these test systems 
we study electronic and structural properties comparing 
the results from the new functional with those obtained 
from "standard" DFT (GGA) and DFT+U calculations. 
Improvements obtained with the novel formulation and 
still remaining issues will be highlighted in each case. 

The paper is organized as follows: in section II the ex- 
tended DFT+U+V energy functional is introduced and 
discussed. In section III we present the results obtained 
from the application of this novel approach to the study 
of bulk NiO, Si and GaAs. Finally, some conclusive re- 
marks are proposed. 



II. THE EXTENDED DFT+U+V FUNCTIONAL 

The LDA+U approach was formulated and developed 
by Anisimov and coworkers in the first part of the 1990s 
[3, Q to improve the accuracy of LDA in describing sys- 
tems characterized by localized, strongly correlated va- 
lence electrons. This scheme is based on a correction to 
the DFT energy functional that can be generally written 
as follows: 

Edft+u — Edft + Eu = Edft + Esub — E^ c . (3) 

In this equation Euub is the part that contains electron- 
electron interactions as modeled in the Hubbard Hamil- 
tonian. Edc is a mean-field approximation to Eu u b and 
models the amount of electronic correlation already con- 
tained in Ed ft- This term is subtracted from the total 
functional to avoid double-counting of the electronic in- 
teractions contained in E^ub- 

Following a commonly used simplified form [28j] of the 
rotationally invariant formulation of DFT+U introduced 
by Liechtenstein et al., [2!| the total corrective functional 
results: 

Eu = E Hub -E dc = J2 \tt [n /CT (l - n /CT )] . (4) 

I,<7 
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In this equation n 1 ' 17 is the "on-site" occupation matrix 
that is defined by the projection of the occupied Kohn- 
Sham orbitals of spin a on the localized (atomic) states 
of orbital quantum number I (typically of d or / kind) 



= ^fkv('<Pkv\ < t>m')(<f>m\'*Pkv)- 



(5) 



In Eq. © m is the magnetic quantum number associated 
with I (—1 < m < I), while /&„ are the occupations of 
the Kohn-Sham orbitals ip% v based on the distribution of 
their energies around the Fermi level. 

The definition of the novel DFT+U+V functional can 
be obtained from a formal generalization of the atomic 
orbital occupation matrix in Eq. ([5]), to allow for the 
possibility that the two atomic wavefunctions involved in 
the definition belong to different atoms: 
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In Eq. © the indexes m and m' run over the angular 
momentum manifolds that are subjected to the Hubbard 
correction on atoms I and J respectively. It is important 
to notice that the occupation matrix defined in Eq. © 
contains informations about all the atoms in the same 
unit cell and the on-site occupations defined in Eq. © 
correspond to its diagonal blocks. 

Using the occupation matrix of Eq. © the energy func- 
tional in Eq. © can be generalized as follows: 
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where, specifically, 
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In Eq. , U 1 and V IJ represent on-site and inter-site 
interaction parameters respectively. The star in the sec- 
ond sum operator reminds that for each atom /, index J 
covers all its neighbors up to a given distance (or belong- 
ing to the same shell) and the trace operator Tr is under- 
stood as follows: Tr (n^n^) - £ , n*£,n#* . 



In Eqs. © and @, 



Ji 



J2„ Tr [n IIa ] is the 



total number of electrons on the atomic states of site I. 

Having been obtained as the extension of a sim- 
plified model (Eq. ©), the inter-site term of the 
DFT+U+V corrective functional (Eq. 0) necessarily 
corresponds to a quite drastic approximation of the 
full off-site many-body part of the electronic Hamil- 
tonian. In fact, terms involving states belonging to 
three or four distinct atoms are all neglected. In 
addition, among the two-site terms only a restricted 
subset is explicitly included in the extended corrective 
Hamiltonian. If we establish the notation: Al/Jf L — 



ijkl 



v Rl 



IK.a JL.tr' r 

n iU rijt - d, 



IL,tr JK.tr' 
trtr'Tl^i LI, 



where V ee corresponds to the kernel of the electronic 
interaction, and <j>l represents the state i of atom /, 
the Coulomb part of the term included in our inter-site 
functional corresponds to A{j/j J (averaged over the 
orbitals of the two sites) . Other contributions describing 
cross charge exchanges between neighbor sites (Aj^ 1 ), 
double (parallel) electron transfers from one site to a 
neighbor one (Af^y), couplings between on-site charge 
and hopping from/to the same site AjJ,^,, and coupled 
intra-site charge exchanges (as, e.g., A^f£) are all 
neglected. The one included in our corrective functional 
(proportional the product of number operators on 
different sites) is thus the only one that has a classical 
Coulomb counterpart and is, presumably, the most 
relevant one. This ansatz, not explicitly tested in this 
work, is indeed consistent with the approximation used 
in the formulation of the model 0, 0] and in the very 
abundant literature where the extended Hubbard model 
is used. 

As evident from Eq. (j4|), the on-site term of the to- 
tal energy introduces a finite energy cost for fractional 
atomic occupations (if \\£ is one eigenvalue of n IIa , Ejj 
is for A„ = or A^ n = 1, positive otherwise); this 
penalty favors a Mott-like ground state with correlated 
electrons localized on atomic orbitals. The effect of the 
inter-site interaction can be easily understood from the 
contribution of the corrective functional in Eq. |JJ to 
the total Kohn-Sham potential. This quantity can be 
computed as the functional derivative of the energy with 
respect to (ipk v )* : 

t r \ i tr \ SE Hub 



= ^ ^ (8mm 
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(10) 



From Eq. (|10p it is evident that the on-site term of 
the potential is attractive for occupied states that are, 
at most, linear combinations of atomic orbitals of the 
same atom (resulting in on-site blocks of the occupation 
matrix, n IIa , dominant on others), whereas the inter- 
site one stabilizes states that are linear combinations of 
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atomic orbitals belonging to different atoms (e.g., molec- 
ular orbitals, that lead to large off-site blocks, n JIrT , of 
the occupation matrix). Thus, a competition sets be- 
tween two opposite tendencies that allows for more gen- 
eral localization regimes and increases the coupling be- 
tween orbitals on different sites. Obviously, the charac- 
ter of the electronic ground state depends critically on 
the relative strength of the on-site (U) and the inter- 
site (V) electronic interactions. The detailed balance be- 
tween these quantities is guaranteed by the possibility 
of computing both parameters simultaneo usly through 
the linear-response approach described in [301 ] . In fact, 
the inter-site interaction parameters correspond to the 
off-diagonal terms of the interaction matrix defined in 
Eq. 19 of Ref. 0. 

It is important to notice that the trace operator in the 
on-site functional guarantees the invariance of the energy 
under rotations of atomic orbitals on the same atomic 
site. However, in the inter-site term, the trace applied 
to the product of generalized occupation matrices, in- 
cluding orbitals from multiple sites, is not sufficient to 
strictly enforce the invariance with respect to rotations 
involving all atomic orbitals in the system at the same 
time (i.e., mixing atomic orbitals that are centered on 
different sites). In fact, to obtain this covariance, the el- 
ements of the electronic interaction matrix should trans- 
form as quadruplets of atomic orbitals and thus have full 
site- and orbital- dependence. At the present level of ap- 
proximation, interactions between more than two atomic 
sites are assumed to be less important than the ones in- 
cluded in Eq. ([7]) and orbital dependence is totally ne- 
glected (the effective interactions included in the present 
formulation can be seen as atomic averages of the fully 
orbital-dependent quantities). This assumption, consis- 
tent with the one operated in the original formulation of 
the extended Hubbard model [g, [7[ , appears reasonable 
for the basis set of atomic orbitals, but may not be gener- 
ally valid if a different basis is used to represent localized 
states. The inclusion of inter-site couplings is, however, 
a significant step towards general covariance with respect 
to the on-site functional. Furthermore, the possibility of 
evaluating both U and V from linear-response theory (30j 
(at least for interactions between couples of atomic sites), 
guarantees a high level of consistency between the atomic 
orbital basis set and the interaction parameters used in 
the functional, and reduces the dependence of the results 
on the specific choice of the localized basis. Site- and 
orbital-dependence of the corrective functional is implic- 
itly included in Wannier-function-based implementations 
of the on-site DFT+U [3l| if the effective interactions are 
computed consistently. However, even in a (e.g., maxi- 
mally localized) Wannier-function-based formulation, no 
guarantee exists that inter-site or multiple-site terms of 
the corrective functional can be neglected (in fact the 
Wannier basis is not generally constrained to be the one 
on which the interaction is site-diagonal). With respect 
to this approach, a formulation based on atomic states 
has the advantage of not requiring the preliminary calcu- 



lation of the localized basis orbitals and offers a physically 
transparent way to select a minimal number of relevant 
interaction parameters to be computed (e.g., based on 
distances between neighbors). 

In the implementation of Eq. (J7J) we have added the 
possibility for the corrective functional to act on two I 
manifolds per atom as, for example, orbitals 3s and 3p in 
Si, or orbitals 4s and 3d in Ni. To the best of our knowl- 
edge, this feature has been implemented only recently 
in a "standard" (on-site) DFT+U functional [32| where, 
however, at variance with our formulation, no interaction 
was established between the two manifolds of the same 
atom. If we call "standard" the higher I states of each 
atom the Hubbard correction acts on, and "background" 
the other I manifold being corrected, for a given couple of 
atomic sites I and J (Eq. 0), we have four interaction 
parameters: standard-standard, standard-background, 
background-standard, and background-background. Al- 
ternatively, the different Hubbard-corrected manifolds on 
each atom can be seen as belonging to different atoms lo- 
cated at the same crystallographic site. The motivation 
for this extension consists in the fact that different man- 
ifolds of atomic states may require to be treated on the 
same theoretical ground in cases where hybridization is 
relevant (as, e.g., for bulk diamond whose bonding struc- 
ture is based on the sp 3 mixing of s and p orbitals). The 
particular choice to have different interaction parame- 
ters acting on different manifolds of atomic states, and 
the possibility to include "cross-manifold" interactions 
(basically "on-site" V parameters between standard and 
background states) eliminates the need for "ad-hoc" con- 
structions of the basis set for cases where localization is 
expected to occur on states different from atomic orbitals 
(see, for example, 33 1, (HI, [HI). 

While the inter-site interaction V can be calculated at 
the same time of the on-site U and with no additional 
cost, the computational workload of a DFT+U+V cal- 
culation depends on the number of neighbors between 
which the inter-site correction is established. For in- 
teractions between nearest neighbors (all contained in a 
3x3x3 supercell centered around the primitive one) 
the computational cost is only marginally bigger than 
for an on-site-only DFT+U; if further shells of neighbors 
are included in the summation of Eq. [7J larger supercells 
around the primitive one may be needed and the search of 
equivalent neighbors may increase significantly the cost 
of the calculation (up to a 18% increase in the cpu time 
for a 5 x 5 x 5 supercell compared to a 3 x 3 x 3). 



III. RESULTS AND DISCUSSIONS 

In this section we present the results obtained with the 
novel functional applied to the study of the selected test 
systems: NiO, Si, and GaAs. While the choice to test a 
numerical approach explicitly designed for strongly corre- 
lated materials on band insulators could appear unusual, 
it is important to notice that, for both classes of mate- 
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rials, the fundamental gap has largely the same nature: 
it is mostly due to the derivative discontinuity of the 
exchange-correlation energy functional (see, for example, 
[iSHS) that is missing from standard DFT approxima- 
tions. Thus a correction, as DFT+U, ideally designed to 
introduce this feature in the energy spectrum of a system, 
should have general validity for semiconductors and in- 
sulators in both families. Analogous considerations have 
led the authors of [38[ to operate a similar selection of 
systems to test their Reduced Density Matrix Functional 
(RDMF) approach. 

The use of the novel functional to study NiO is moti- 
vated by the need to validate it on systems (as Mott or 
charge-transfer insulators) for which the on-site DFT+U 
already gives reasonably good results. Also, the charge- 
transfer character of the band gap of this material (i.e., 
the hole in the O ligand field observed in photoemission 
experiments) |30( ] suggest a more significant hybridiza- 
tion between d and p states than in other oxides. Si and 
GaAs are, instead, typical band insulators, with cova- 
lent and ionic character, respectively. The application 
to these materials is a very important one for the novel 
functional introduced in this work as the electronic struc- 
ture of these systems is based on the formation of bonds 
and thus on the hybridization between states of neighbor 
atoms. 

For all these systems comparison will be made with the 
results obtained with DFT and with the "standard" , on- 
site DFT+U. All the calculations were performed using 
the plane-waves pseudo-potential "pwscf ' code contained 
in the Quantum ESPRESSO (QE) package [H EH, 
where we implemented the "+U+V" correction starting 
from an existing on-site DFT+U functional [4l[. Ex- 
cept in one case (detailed in the paper), all the pseu- 
dopotentials were obtained from the QE pscudopotcntials 
database [io| . 

It is important to remark that for all the systems 
treated in this work U and V were obtained using the 
linear response approach described in [301 ] - However, at 
variance with the calculations presented in [30l ]. the re- 
sponse matrices were not constrained to give a neutral 
total response (in terms of on-site occupations) when the 
perturbation is applied on every atom. While this has 
no effects on the final results if big enough supercells 
are used in the calculations, we believe it to be a bet- 
ter strategy to compute the effective interactions espe- 
cially if the response of all atomic orbitals (standard plus 
background) is explicitly considered. A further difference 
with the results presented in [30( consists in the use of 
orthogonalized atomic orbitals to construct occupation 
matrices. While not necessary, this choice guarantees 
that the atomic occupation matrices satisfy more strin- 
gent sum rules (their trace is closer to the total number 
of electrons on the atomic states). Some small differences 
are to be expected with the results obtained for NiO in 
[3fj| ; however, we believe that the consistent evaluation of 
the effective electronic interactions (especially if obtained 
from a DFT+U ground state) reduces these differences 



to a minimal value. 



A. NiO 

As other transition-metal oxides with the same stoi- 
chiomctry, nickel oxide has a cubic rock-salt crystalline 
structure. However, below a Neel temperature of 523 K, 
the magnetic moments of Ni atoms arrange in an anti- 
ferromagnetic (AF) order (usually referred to as AFII) 
where ferromagnetic (Iff) Ni planes alternate with op- 
posite magnetization. As a consequence of the AF mag- 
netic order the crystal acquires a rhombohedral symme- 
try. Transition metal oxides (TMOs) have represented 
a significant challenge for theorists since their insulat- 
ing character can not be explained satisfactorily using 
band theory due to the intrinsic many-body origin of 
their band gap. Although within Kohn-Sham theory 
the energy spectrum has no precise physical meaning 
and is not bound to reproduce the band gap of the sys- 
tem, once the finite discontinuity of the (exact) exchange- 
correlation potential is added to the HOMO-LUMO en- 
ergy difference the correct band should be exactly repro- 
duced [H, [37j|- Failing to properly incorporate many- 
body effects, most of approximate exchange-correlation 
functional produce no discontinuity in the corresponding 
potentials, thus resulting in quite poor estimates of band 
gaps or in their total suppression. In fact, many TMOs 
are predicted to be metal in striking contrast with the 
observed insulating character. In some of these systems 
(e.g., FeO [30() lack of proper account of many-body as- 
pects of the electronic structure also results in serious in- 
accuracies in the description of structural and magnetic 
properties. Similar problems can be expected for NiO 
even if a gap (still significantly smaller than the exper- 
imental one) fortuitously opens in its Kohn-Sham spec- 
trum as the result of the balance between the exchange 
and the crystal-field splittings among the d states of Ni 
(nominally occupied by eight electrons). Photoemission 
experiments on NiO have measured a band gap of about 
4.3 eV (3.1 eV at the minimum intensity) [39( and have 
explored several features of the spectrum asses sing , in 
particular, its charge-transfer character [3!| HH, EJ ^ n 
fact, the excitation of one electron across the band gap 
of the quasi-particle spectrum corresponds to its trans- 
fer from the p states of a ligand oxygen atom to the 
d states of a Ni atom. As a consequence, the top of 
the valence band is dominated by the p states of oxygen 
while the bottom of the conduction one largely consists 
of Ni d states. Computational studies on this material 
have been quite successful in reproducing these features 
of its spectrum. DFT+U, whose corrective potential is 
designed to introduce a finite energy difference between 
occupied and unoccupied states, can be expected to pro- 
duce a more accurate estimate of the fundamental gap 
of a system (in principle larger than the one appearing 
in photoemission experiments). In fact, it has been used 
quite successfully to study NiO and has produced a band 
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gap of about 3.0 to 3.5 eV (the precise value varies with 
the U used in different works) and quite accurate esti- 
mates for both the magnetic moments and the equilib- 
rium lattice parameter [44], EH, Efil]- For other details 
of the density of states (as, e.g., the spectral weight of 
O p states on top of the valence band), the agreement 
was not unanimous. DFT+U has also been employed re- 
cently to compute the fc-edge XAS spectrum of NiO using 
a novel, parameter-free computational approach [47| that 
has produced results consistent with experimental data. 
In this articles authors highlight the importance of non- 
local excitation of d states and, specifically, those involv- 
ing second-nearest neighbor Ni atoms. The same intent 
of improving the description of spectroscopic properties 
has been pursued with the GW approximation based on 
a DFT+U ground state [HI, [l3j] ; this approach has pro- 
vided a better estimate of the energy gap compared to 
DFT+U even though other details of the density of states 
were almost unchanged 12). Hartree-Fock (HF) |48[ and 
hybrid functional (e.g. B3LYP) 0, El have also been 
used to study NiO. While pure HF overestimates the size 
of the band gap and fails to reproduce its charge-transfer 
nature, hybrid functionals are more accurate and result in 
a better estimate of the properties of this material. Most 
recently DFT+DMFT calculations of NiO have explored 
the effects of dynamical correlations and have produced 
a band gap in the energy spectrum in excellent agree- 
ment with photoemission experiments, even though its 
charge-transfer character was not always well reproduced 
[HoL l5l| . A modification to this method within the itera- 
tive perturbation theory has been recently used to fix this 
particular aspect of the calculated spectrum of NiO [H[ . 
In this paper, the authors also present a detailed discus- 
sion of the features of the DOS remarking the role of 
dynamical correlations (e.g., in eliminating the excessive 
spectral weight of mixed Ni e g and O p states at the bot- 
tom of the valence band, normally observed in DFT+U 
and GW), and the importance that inter-spin and inter- 
atomic components of the self-energy would have in fix- 
ing some still remaining discrepancies in the computed 
spectrum compared to the experimental one (e.g., the 
position of the occupied satellite). This latter hypothe- 
sis will be discussed with the results of the present work 
where inter-atomic couplings are explicitly taken into ac- 
count thanks to the extended functional introduced in 
the previous section. 

All calculations presented in this paper were performed 
in the AF configuration of the system. We used ul- 
trasoft [53| GGA pseudopotentials constructed with the 
PBE parametrization [5J| (Ni.pbe-nd-rrkjus.UPF and 
O . pbe-rrk j us . UPF of the Quantum-Espresso repository 
[40]) and based on the following valence configurations: 
3rf 8 4s 2 4p° for Ni and 2s 2 2p 4 for O. Energy cut-offs of 
40 and 400 Ry were used for wavefunctions and charge 
density respectively. A 4 x 4 x 4 Monkhorst-Pack fc-point 
grid [55J was used to sample the Brillouin zone of this sys- 
tem. Due to its magnetization, NiO was formally treated 
as a metal with a 0.01 Ry gaussian smearing of its Fermi 
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FIG. 1: (Color online) The variation of the on-site and inter- 
site electronic couplings as functions of the cubic lattice pa- 
rameter of NiO. The interaction parameters were obtained 
from a GGA ground state (top graph) or from a GGA+U+V 
ground state (bottom graph) according to a self-consistent 
procedure (see text). Black lines represent on-site (Ni only) 
interactions, red ones describe the couplings between nearest 
neighbors (Ni-O), blue ones between next-nearest neighbors 
(Ni-Ni). Solid lines are for standard-standard interactions, 
dashed for standard-background, dot-dashed for background- 
background. 



distribution function that conveniently placed the (fake) 
Fermi level in the middle of the energy band gap. 

The electronic interaction parameters (U and V) used 
in this work were re-computed for each considered lattice 
spacing (using the linear-response approach of ref. [30l| ) . 
Fig. [1] shows the dependence of some of these parameters 
on the size of the unit cell. It is important to remark that 
these parameters are not directly comparable with those 
evaluated in [3fl | for two reasons: first, in the present 
calculations we have explicitly included the response of 
the "background" states, which could not be done in the 
previous work; second, at variance with what was done in 
Ref. [3(| to define the on-site occupation matrix we have 
used atomic orbitals that were preliminarily orthogonal- 
ized. However, since Us and Us are consistently com- 
puted and used within the same approximation (i.e., us- 
ing the same definition of the occupations) in both works 
the most significant differences are expected to descend 
from using a U-augmented functional. 
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As it can be observed from Fig. Q] (top graph) the 
on-site U is the parameter that changes the most with 
the size of the unit cell. Other electronic couplings are 
much less affected and the dependence of the interaction 
between second-nearest neighbors (Ni-Ni) on the lattice 
parameter is almost negligible. While the parameters 
shown in the top graph of figure [T] were obtained from 
the GGA ground state of the system, those represented 
in the bottom graph were computed "self-consistently" 
from a GGA+U+V ground state using a method equiv- 
alent to that of ref. (5(| but based on a more efficient 
algorithm that will be described elsewhere. Although 
all the electronic interactions are obtained at the same 
time (including those for "background" states and their 
coupling with "standard" ones) the interaction param- 
eters between the p states and between p and s states 
of oxygen (all Uo and Vb-o) were excluded from our 
DFT+U+V calculations. While a rigorous motivation 
for this choice is presently missing, we believe that the 
corrective action of the "+U+V" functional is mostly 
needed for open manifolds around the Fermi level (or the 
top of the valence band) while p and s states of oxygen 
are almost completely full. Besides the on-site C/jvi> first 
(Vni-o) and second (Vm-m) nearest neighbor inter-site 
couplings (between "standard" and "background" states) 
were all included in our calculations; the inclusion of the 
interactions between further neighbors was found to pro- 
duce no relevant effect. For all the approaches compared 
in this work (GGA, GGA+U and GGA+U+V) the op- 
timization of the rock-salt cubic structure (in the AF 
ground state) was performed using the cell-size depen- 
dent electronic couplings plotted in Fig. [TJ fitting the 
the dependence of the total energy on its volume on 
a Murnaghan equation of state. Using size-dependent 
electronic couplings has been recently shown to be fun- 
damental for quantitative descriptions of structural and 
electronic properties of materials and for accurate evalu- 
ations of transformation pressures between, e.g., different 
spin or structural phases [57], [58| . 

Table [J compares the equilibrium lattice parameters, 
the bulk moduli, and the energy band gap obtained in 
the three different approaches with results from experi- 
ments. In each case the band gap of the material was 
measured for the equilibrium lattice parameter reported 
in the same table. A comparison between the DOS ob- 
tained within the different approaches is made in Fig [21 

Our GGA calculations result in an equilibrium lattice 
spacing of 7.93 bohrs which is quite consistent (+0.5%) 
with the experimental value of 7.89 bohrs and provides 
a good estimate of the bulk modulus (187.7 GPa). As 
expected from other GGA results, the band gap (~ 0.6 
eV, Fig. [21 top graph) is smaller than the observed one 
(3.0 - 4.3 eV) and its spectroscopic nature is also not con- 
sistent with experiments as the top of the valence band 
is dominated by the d states of Ni instead of O p. It is 
important to notice that the peak at 5 eV below the edge 
of the valence band contains contributions from both O 
and Ni and, as specified in other works in literature (see, 




Energy (eV) 



5 I — ■ — i — ■ — i — ■ — i — ■ — i — ■ — i — ■ — i — ■ — r 



CD 




Energy (eV) 



4 I — ■ — i — ■ — i — ■ — i — ■ — i — ■ — i — ■ — i — ■ — r 

"55 

o 




Energy (eV) 



4 I — ■ — i — ■ — i — ■ — i — ■ — i — ■ — i — ■ — i — i — r 

O 

:> 3 - 

CD 




Energy (eV) 



FIG. 2: (Color online) The Density of states of NiO obtained 
with different approximations: a) GGA; b) GGA+U (U on Ni 
d states only); c) GGA+U including on-site interactions on 
Ni d (Uni) and O p (Uo) states; d) GGA+U+V including on- 
site (Ujvi and Uo) and inter-site (Vni-Ni, Vo-o, and Vni-o) 
interactions. Only "standard" states corrected. In all graph 
the black and red lines represent majority and minority spin 
d states of Ni, the blue line the p states of O. 
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TABLE I: The equilibrium lattice parameter, (a, in Bohr 
atomic radii), the bulk modulus (B, in GPa), and the band 
gap (Eg, in eV) of NiO obtained with different computa- 
tional approaches: GGA, "traditional" GGA+U (with U only 
on the d states of Ni), GGA+U+V and a "self-consistent" 
GGA+U+V with the interaction parameters computed from 
a GGA+U+V ground state (see text). Comparison is made 
with the experimental results on all the computed quantities. 





a 


B 


E g 


GGA 


7.93 


188 


0.6 


GGA+U 


8.069 


181 


3.2 


GGA+U+V 


8.031 


189 


3.6 


GGA+U+V SC 


7.99 


197 


3.2 


Exp 


7.89" 


166-208 6 


3.1-4.3 C 



"Reference M 
'Reference 6(1 
Reference |39 t E^l 



e.g., [il]), is the signature of the hybridization between 
d eg and p states. This is a spurious feature, common 
also to DFT+U calculations, as the lowest part of the 
valence state should be dominated by ti g states instead 

GGA+U (U on Ni d states only) leads to a significant 
improvement in the size of the band gap (3.2 eV) com- 
pared to GGA and O p states are dominant on the top 
of the valence band in agreement with experimental re- 
sults (Fig. [5] second graph from the top). The equilibrium 
lattice parameter obtained with GGA+U (8.069 bohrs) 
is larger than the one obtained with GGA, deteriorat- 
ing the agreement with experiments (table HJ . Also, the 
crystal results softer with a bulk modulus of about 181 
GPa (still in the range of experimental results). While 
the almost exclusive presence of p bands on top of the 
valence band seems to overestimate the spectral weight 
of these states compared with the results of other cal- 
culations [3 |4f| [40|, the dominant p character of the 
peak at -2 eV, and the appearance of a peak at -4.5 eV, 
mainly due to minority-spin t2a Ni states, have not been 
highlighted in experiments [43j. However, the position 
of a strong peak in the DOS at about -7 eV from the 
valence band edge appears consistent with experimental 
results [39] although further analysis would be needed to 
assess its character. The clear-cut separation between p 
states and lower energy d states is, instead, unusual for 
DFT+U calculations. This latter aspect may be due to 
the different value of the U parameter used with respect 
to other calculations of this kind. 

Using GGA+U+V, the equilibrium lattice parameter 
(Table fl} is corrected towards the experimental value 
(8.03 bohrs) while the bulk modulus results equal to 189 
GPa and a band gap of about 3.6 eV opens between O p 
and Ni d states correctly placed at the top of the valence 
band and the bottom of the conduction one respectively. 
The band structure of the system shows some differences 
from the one obtained with GGA+U. With respect to 



the GGA+U results, the minority spin d states of Ni 
have moved upwards in energy and show a larger overlap 
with the p states of O that suggests a larger degree of 
hybridization. As a result, the peak at -2.5 eV acquires a 
markedly mixed d and p character, and the central peak, 
mainly due to Ni minority-spin d states, moves to slightly 
higher energy at -4 eV from the gap. The occupied satel- 
lite peak, that also appears at slightly higher energy (« 
-6.5 eV) than in the GGA+U DOS, is still dominated by 
Ni majority-spin d states. 

Using the "self-consistent" GGA+U+V (i.e., with U 
and V computed from a GGA+U+V ground state) we 
obtained an equilibrium structure in better agreement 
with experimental results (a = 7.99 bohrs), a bulk mod- 
ulus of about 197 GPa (still within the range of available 
experimental data) and a band gap of 3.2 eV (at the ba- 
sis of the valence and conduction peaks) as evident from 
the bottom graph of Fig. O The electronic structure 
in this case is quite similar to that obtained within the 
non self-consistent GGA+U+V approach (Fig. [2j third 
graph from top). However an even larger overlap can be 
observed between Ni d and O p states in the upper part 
of the spectrum with a peak below the edge of the va- 
lence band at w -2.2 eV that is dominated by Ni states in 
agreement with experiments 39] . The occupied satellite 
peak, at about -6 eV, is still dominated by Ni d states 
although the p state component seems slightly more sig- 
nificant in this case. It is important to notice that the 
valence DOS has now acquired three dominant features 
that correspond to the O peak at the top of the band, 
a Ni-dominated peak at about -2.2 eV and a strong d 
peak at -6 eV. These characteristics are consistent with 
the results from photoemission experiments although the 
relative spectral weight of these peaks are not in quan- 
titative agreement with observations. Furthermore we 
don't observe the downshift of the occupied satellite (ac- 
tually moved to slightly higher energies) that in Ref. [52] 
is expected to descend from the explicit inclusion of inter- 
site electronic interactions in the corrective functional. In 
our opinion these discrepancies are due to the fact that, 
at the present level of approximation, electronic interac- 
tions parameters are only site-dependent. While dynam- 
ical correlation may play an important role, as demon- 
strated by DFT+DMFT calculations [U EH S3 , we be- 
lieve that the most significant refinement needed to re- 
produce all the details of the spectrum of the material 
consists in the formulation of an extended model with 
orbital-dependent interaction parameters. Specifically, a 
distinction should be made between the (inter-site) inter- 
actions of O p states with Ni ti g and e g states. 

In summary, while not completely satisfactory, 
GGA+U+V offers a quite significant improvement in the 
description of NiO compared to GGA+U and we believe 
that the more realistic representation of the electronic 
properties (as well as structural ones) can be a better 
starting point than GGA+U for numerical approaches 
(as GW or DFT+DMFT) aimed at predicting the exci- 
tation spectrum of the material or at studying the effects 
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of dynamical correlations. 



B. Si and GaAs 

For Si and GaAs, DFT calculations based on LDA 
and GGA approximate functionals provide results for 
the structural properties (e.g., the lattice parameter and 
the bulk modulus and its first derivative) [6ll [6^, [fiH 
and the vibrational spectrum [64[ in good to excellent 
agreement with experiments. GGA normally results in 
slightly longer equilibrium lattice constants and softer 
bulk moduli than the measured ones [13, HH, [H| . How- 
ever, the computed band gap of these materials systemat- 
ically underestimates the experimental value due to inac- 
curacies inherent to approximate DFT energy functionals 
and, in particular, to the lack of finite discontinuities in 
the exchange-correlation potential [H, [37} . Many correc- 
tive approaches have been developed to overcome these 
difficulties of standard DFT approximations and quite 
good agreement with the experimental spectrum could 
be obtained from numerical approaches correcting the 
exchan ge c orrelation energy model, as in exact-exchange 
( EXX) [661 |67| or hybrid functional [68| calculations, and 
from methods specifically designed to reproduce the elec- 
tronic excitation spectrum, as the GW approach based 
on an LDA [H, E3] or a EXX [ED] ground state. In 
the present work, we investigate the performance of the 
GGA+U+V functional on Si and GaAs and compare its 
results with those from standard GGA and GGA+U. In 
these semiconductors, the hybridization between s and 
p orbitals and the formation of covalent bonds between 
neighbor atoms are key ingredients to understand their 
structural and electronic properties. Thus, the correc- 
tive functional introduced in the present paper, includ- 
ing inter-site interactions and acting on two angular mo- 
menta manifolds per atom, contains all the necessary in- 
gredients to achieve significant improvement with respect 
to standard DFT approximations in the description of 
these materials. 

As mentioned above, for both Si and GaAs we have 
used GGA exchange-correlation functionals constructed 
with the PBE parametrization [54J. In particular for Si 
we employed a norm-conserving pseudopotential (Si.pbe- 
rrkj.UPF of the QE repository) that was built on the 
3s 2 3p 2 3d° valence configuration. This potential required 
an energy cut-off of 40 Ry for the plane-wave expansion 
of the Kohn-Sham wavefunctions. For GaAs we used, in- 
stead, ultrasoft pseudopotentials [53| that required cut- 
offs of 40 and 400 Ry for the expansion of the electronic 
wavefunctions and charge density respectively. For As 
we chose a pseudopotential constructed on the 4s 2 4p 3 
configuration (As.pbe-n-van.UPF of the QE library). As 
reference [72j points out, predictions for gallium com- 
pounds can be very sensitive to the inclusion of 3c? elec- 
trons in the valence manifold for Ga pseudopotentials. 
Thus, two alternatives were considered for Ga with 3d 
electrons in valence (3d 10 4s 2 4p 1 - Ga.pbe-nsp-van.UPF 
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FIG. 3: Intra-site (U) and inter-site (V) interaction param- 
eters as functions of lattice parameter (a) for GaAs with 
3d electrons as valence electrons, a) Intra-site interactions. 
Squares correspond to Ga atom and circles to As atom. 
Filled symbols represent U 3S (standard-standard), open sym- 
bols represent U ab (standard-background) and symbols with 
crosses represent U bb (background-background), b) Inter-site 
interactions. Arrows indicate the equilibrium lattice parame- 
ter. 



TABLE II: Interaction parameters U and V (eV) for Si and 
GaAs 6 . Inter-site terms are for first-neighbors and the listed 
values are for the equilibrium lattice parameters found with 
GGA+U+V. Indexes s and 6 stand for "standard" and "back- 
ground" orbitals respectively as discussed in the text. 
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1.68 
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'Ga 3d electrons as valence electrons. 



from the QE repository), or frozen in the core [As Ap ). 
We will refer to these two different situations as GaAs(v) 
and GaAs(c) respectively. The Ga pseudopotential for 
GaAs(c) was constructed using D. Vanderbilt's code and 
library [73|]. For both materials (Si and GaAs), the Bril- 
louin zone was sampled with a 8 x 8 x 8 Monkhorst-Pack 
[55l | k-point grid. 

In order to determine the equilibrium lattice parameter 
and the bulk modulus, the dependence of the ground- 
state energy on volume was fitted by the Murnaghan 
equation of state. As in the case of NiO, DFT+U and 
DFT+U+V total energy calculations were performed us- 
ing interaction parameters consistently recomputed for 
each value of the lattice parameter. Fig. [3] illustrates the 
dependence of Us and Vs on the size of the unit cell in 
the case of GaAs(v). 

In Table UH the on-site and inter-site interactions com- 
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TABLE III: Comparative results for lattice parameter (a, in 
A), bulk modulus (B, in GPa) and energy gap (E g , in eV). 

Si GaAs(v) a GaAs(c)'' 

a B E g a B E a a B E g 

GGA 5.479, 83.0, 0.64 5.774, 58.4, 0.19 5.578, 65.7, 1.25 
+U 5.363, 93.9, 0.39 5.736, 52.6, 0.00 5.616, 62.7, 0.81 
+U+V 5.370, 102.5, 1.36 5.654, 67.7, 0.90 5.535, 76.5, 1.97 

Exp. c 5.431, 98.0, 1.12 5.653, 75.3, 1.42 5.653, 75.3, 1.42 

a Ga 3d electrons as valence electrons. 
b Ga 3d electrons in the core. 

c at 300 K, from Ref. Q- At K, the energy gaps of Si and GaAs 
are estimated to be 1.17 eV and 1.52 eV respectively [74)1 . 



puted at the equilibrium lattice parameter are grouped 
to ease the comparison. It is important to notice that 
the values of either the Us or the Vs are very similar ir- 
respective of the orbital manifold the parameters acts on 
(p or s states). This result, more evident for GaAs and 
for the inter-site couplings, is due to the hybridization 
of p and s orbitals and corroborates the need of treating 
both orbitals at the same level. 

In Table IIII1 the equilibrium lattice parameter, the 
bulk modulus and the band energy gap obtained from 
GGA, GGA+U and GGA+U+V calculations on Si and 
GaAs can be directly compared with experimental mea- 
surements of the same quantities (we refer to the data col- 
lected in the web-database, Ref. [ZJ])- As expected from 
abundant literature, GGA overestimates the equilibrium 
lattice parameter with respect to the experimental value 
(except for the case of GaAs(c)), while the bulk modulus 
and the band gap are underestimated in all cases. On- 
site GGA+U predicts the equilibrium lattice parameter 
in better agreement with the experimental value (over- 
corrected for Si); however, the bulk modulus is improved 
with respect to the GGA value only in the case of Si. In 
all three cases, however, the energy band gap is lowered 
compared to GGA, further worsening the agreement with 
experiments. While this result may appear strange, it is 
not unexpected. In fact the on-site corrective functional 
suppresses the hybridization between states on neigh- 
bor atoms that is largely responsible for the band gap 
in semiconductors between valence (bonding) and con- 
duction (anti-bonding) states. The use of the inter-site 
correction, in spite of the fact that Vs are smaller than 
half the on-site Us, results in a systematic improvement 
for the evaluation of all these quantities. In fact, en- 
couraging the occupations of hybrid states, the inter-site 
interactions not only enlarge the splittings between pop- 
ulated and empty orbitals (which increases the size of 
the band gap), but also make bonds shorter (so that hy- 
bridization is enhanced) and stronger, thus tuning both 
the equilibrium lattice parameter and the bulk modulus 
of these materials to values closer to experimental results. 
For GaAs the calculations with the Ga 3d states in va- 
lence are more accurate than the ones with these atomic 



states frozen in the core. In fact, for GaAs(v) we obtain 
an equilibrium lattice parameter in excellent agreement 
with the measured one, while bulk modulus and band gap 
slightly underestimate the experimental results, yet im- 
proving the estimates obtained with GGA and GGA+U. 
For GaAs(c), instead, while the bulk modulus is in very 
good agreement with experiments, the equilibrium lattice 
spacing is underestimated and the band gap significantly 
overestimated worsening, for both quantities, the accu- 
racy of the results obtained with GGA and GGA+U. Our 
results thus confirm that Ga 3d should be treated as va- 
lence states; in the present study, however, they are not 
directly subject to the Hubbard functional. 

Despite the overall improvement obtained with 
GGA+U+V, some discrepancies with the experimental 
results still persists regarding, especially, the equilib- 
rium lattice constant and the bulk modulus. However, it 
should be kept in mind that computational results pre- 
sented in this section would be directly comparable to 
K measurements for which slightly shorter lattice param- 
eters and slightly higher bulk moduli are predicted [74| 
(see footnotes of Table IIIII for details) . 

Fig. g] displays the GGA and GGA+U+V band struc- 
tures and density of states for Si and GaAs. The energies 
were shifted so that the top of the valence band corre- 
sponds to the zero of the energy in all cases. The cor- 
rection introduced by the GGA+U+V (with respect the 
GGA band structure) to the conduction band consists in 
an almost rigid upward shift in energy. In the valence 
manifold, instead, the correction to the energy levels ac- 
quires a slight k— dependence and has a more pronounced 
effect on the lowest energy level, which is shifted down- 
ward. This latter effect results in an increase of the total 
bandwidth of these systems (by 0.6 eV for GaAs and 
by 0.8 eV for Si) compared to GGA, and accounts for 
perfect agreement with experimental results (measured 
bandwidth is 12.5 eV for Si and 13.2 eV for GaAs \&^). 

In summary, the GGA+U+V shows significant im- 
provement over GGA and GGA+U in the description 
of structural and electronic properties of Si and GaAs. 
Although the dominant effect arises from the interaction 
between nearest neighbors sites, for both Si and GaAs 
the inclusion of the coupling between further neighbors 
could bring a refinement of the results presented in this 
section. Also, in the case of GaAs with Ga 3d electrons 
in valence, more accurate results could be obtained ac- 
counting for the hybridization of these states with 4s and 
A.p orbitals (probably the overlap with the states of As 
is quite limited) and including the corresponding inter- 
action parameter in the "+U+V" corrective functional. 
Since the GGA and the GGA+U+V ground states are 
qualitatively similar, we expect a "self-consistent" evalu- 
ation of U and V to have a minor effect on the presented 
results both for Si, and GaAs. 
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IV. CONCLUSIONS 
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FIG. 4: Band structure and density of states of Si and GaAs 
(core case). Continuous lines represent GGA+U+V results 
and dashed lines represent standard GGA results. All ener- 
gies were shifted so that the top of valence bands are at zero 
energy. 



In this paper we have introduced a useful general- 
ization (named DFT+U+V) of the popular DFT+U 
method that, modeled on the extended Hubbard Hamil- 
tonian, includes on-site (U) and inter-site (V) electronic 
interactions and is able to correct up to two angular mo- 
mentum manifolds per atom. The competition between 
the two kinds of electronic couplings avoids the over- 
stabilization of occupied atomic orbitals, often affecting 
the on-site DFT+U, and allows for the description of 
more general ground states where electrons "localize" on 
hybridized (e.g., molecular) orbitals. The flexibility in 
the representation of occupied states is further enhanced 
by the larger number of orbital manifolds that are simul- 
taneously subject to the action of the corrective func- 
tional. Numerical accuracy is guaranteed by the linear- 
response calculation of U and V that allows to evaluate 
the relative strength of the two couplings and, thus, to 
precisely determine the degree of electronic localization 
of the ground state. 

The effectiveness of the method is demonstrated in this 
paper by its successful application to quite diverse test 
systems as NiO (Mott/charge-transfer insulator), Si and 
GaAs (band semiconductors). For both classes of ma- 
terials the use of the DFT+U+V functional results in 
a significant improvement (over approximate DFT and 
DFT+U) in the agreement of structural (equilibrium lat- 
tice parameter and bulk modulus) and electronic (overlap 
of p and d states, band gap) properties with available ex- 
perimental results. For Si and GaAs this improvement 
is more consistent, probably due to the higher level of 
structural isotropy of these materials that fits better with 
a corrective functional with orbital independent interac- 
tions. 

The importance of these results is not limited to the 
more accurate description the selected test systems re- 
ceived. In fact, they demonstrate that the novel compu- 
tational approach provides a unified theoretical frame- 
work able to treat systems as divers as Mott, charge- 
transfer, and band insulators with equal accuracy and 
comparable computational effort. Furthermore, the use 
of the extended Hamiltonian provides a new route to im- 
prove the corrective functional in a systematic way as it 
allows to add interactions between further shells of neigh- 
bors independently and, thus, to easily assess their indi- 
vidual contributions to the final results. This procedure 
also leads to the possibility to set calculations of minimal 
cost (with a minimal set of inter-site interaction terms). 

Potential applications of this novel corrective func- 
tional are quite broad. Besides high T c superconduc- 
tors (for which the nearest-neighbor electronic couplings 
has sometimes been demonstrated to play a very impor- 
tant role), the remarkable success obtained with strongly 
localized (correlated) and strongly hybridized systems, 
suggests that intermediate situations where "correlated" 
and "non-correlated" orbitals show significant overlap 
(as, e.g., semiconductors doped with magnetic impuri- 
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ties, metallic active centers of molecular complexes, etc) 
and phenomena where a switch between different local- 
ization regimes is observed (e.g., in bond breaking and 
formation events) are likely to be described quite accu- 
rately by the novel DFT+U+V scheme. This is espe- 
cially true for isotropic systems that better fit the simple 
"+U+V" functional based on orbital-independent inter- 
actions; more anisotropic materials may require, instead, 
extra attention. The development of a more flexible cor- 
rective Hamiltonian based on orbital-dependent effective 
interactions, needed in these cases, is planned for future 
investigations. 
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